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INTRODUCTION 


PURPOSE . 

The  purpose  of  this  project  was  to  develop  and  build  a  prototype  weather 
contouring  device  to  aid  in  the  preparation  of  a  procurement  specification. 

The  device  was  to  be  retrofitted  to  the  Airport  Surveillance  Radar  (ASR-8)  and 
be  usable  on  any  conventional  plan  position  indicator  (PPI).  It  was  to  pro¬ 
vide  two  specific  contour  types  chosen  by  the  Air  Traffic  Service  personnel. 

Each  type  is  capable  of  displaying  two  of  a  possible  six  predetermined  levels 
of  weather  intensity. 

The  contours  produced  by  this  device  were  to  be  related  to  units  of  reflectivity 
(dBZ)  rather  than  simply  received  power  or  some  equally  arbitrary  units.  It 
was  assumed  that  the  input  video  would  be  log  normal,  and  no  attempt  was  made 
to  solve  the  moving  target  indicator  (MTI)  velocity  response  or  MTI  limiting 
problems,  even  though  MTI  must  be  used  in  a  practical  system. 

BACKGROUND . 


The  advent  of  circular  polarization  (CP),  logarithmic  fast  time  constant  (LOG 
FTC),  and  other  techniques  for  eliminating  weather  from  the  air  traffic  con¬ 
trollers'  display  has  made  it  desirable  to  provide  weather  information  to  the 
operator  in  other  than  the  original  raw  radar  form,  specifically  to  provide 
the  air  traffic  controller  with  continuous  weather  information  superimposed  on 
the  air  traffic  control  (ATC)  display  in  a  form  that  does  not  interfere  with 
the  detection  and  control  of  aircraft  targets.  This  display  of  weather  on  the 
controller  display  has  generally  been  called  weather  contouring. 

Weather  contouring  originally  referred  to  the  generation  of  a  thin  isoreflec¬ 
tivity  line  along  the  weather  system  perimeter.  Later,  the  term  was  expanded 
to  include  any  processed  display  of  weather  reflectivity.  Weather  contours 
were  lines  generated  by  displaying  values  of  reflectivity  between  Z  and  Z  +^Z, 
where  Z  was  the  threshold  level  and  \Z  was  an  arbitrarily  small  reflectivity 
increment.  For  small and  a  smooth  weather  function,  this  produced  a  smooth, 
continuous  isoreflectivity  line.  In  reality,  weather  returns  do  not  produce 
a  well-defined,  smooth  function  and  so  do  not  produce  smooth  or  continuous 
contours.  An  increase  in  the  value  of  ^Z  coupled  with  range  smoothing  of  the 
weather  video  signal  causes  the  generated  lines  to  be  continous  and  also  some¬ 
what  wide  and  therefore  tends  to  mask  out  the  discontinuities  of  the  contour. 
This  is  not  without  cost,  since  weather  near  the  threshold  which  does  not  change 
abruptly  will  produce  excessively  wide  lines  adversely  affecting  resolution  and 
cluttering  the  display. 

Elimination  of  the  varying  line  width  could  be  accomplished  with  the  intro¬ 
duction  of  a  single-pulse  generating  circuit  (one-shot)  after  the  threshold 
operation.  This  generates  lines  with  constant  radial  dimension.  However, 
many  of  the  lines  which  are  needed  to  outline  weather  are  in  the  radial  direc¬ 
tion  and  require  radial  dimensions  much  larger  than  the  single-pulse  width. 

The  single-pulse  circuit  produces  discontinuities  in  this  instance. 
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Also,  the  early  weather  contouring  systems  produced  misregistration  in  range 
and  azimuth  caused  by  the  time  delay  introduced  by  the  contouring  process. 

In  response  to  a  request  by  the  Federal  Aviation  Administration  (FAA),  Air 
Traffic  Service  (ATS)  for  improved  weather  detection  and  display  capabilities 
for  the  ATC  radar  system,  the  device  described  in  this  report  was  developed  by 
the  National  Aviation  Facilities  Experimental  Center  (NAFEC)  to  meet  the  need 
for  improved  display  of  weather  information  and  to  solve  the  above  problems. 


DESCRIPTION  OF  EQUIPMENT  REQUIREMENT 


DISCUSSION. 

A  weather  cell  is  composed  of  a  large  number  of  small  targets;  i.e.,  raindrops, 
which  move  slightly  between  radar  samples.  When  the  raindrops  are  positioned 
such  that  the  returned  signals  add  together  mostly  in  phase  with  one  another, 
the  received  signal  is  large.  When  the  returned  signals  add  together  mostly 
out  of  phase,  the  received  signal  is  small.  It  is  this  variation  of  the 
received  signal  that  causes  the  principal  weather  display  problem.  Note  that 
any  given  measurment  of  the  amplitude  of  a  weather  cell  can  vary  from  the  mean 
by  a  considerable  amount;  in  fact,  variations  of  greater  than  5  decibels  (dB) 
can  be  expected  about  half  of  the  time. 

To  make  a  better  estimate  of  the  weather  amplitude,  several  samples  must  be 
averaged  together.  The  greater  the  number  of  samples  averaged,  the  better  the 
estimate.  This  averaging  can  be  accomplished  in  three  ways:  (1)  averaging 
consecutive  pulses  in  range,  (2)  averaging  adjacent  range  segments  in  azimuth, 
and  (3)  averaging  over  a  given  area  for  a  number  of  consecutive  antenna  scans. 
Range  and  scan  averaging  increase  the  number  of  effective  independent  samples 
by  the  actual  number  of  samples  averaged,  while,  because  of  the  signal  correla¬ 
tion  in  azimuth,  azimuth  averaging  increases  the  number  of  effective  samples 
by  a  factor  less  than  the  number  of  actual  averaged  samples.  Both  range 
and  azimuth  averaging  reduce  the  spatial  resolution  of  the  weather  display. 
Typically,  the  accuracy  of  the  weather  display  is  determined  by  the  required 
resolution  and  the  resolution  inherent  in  the  radar  system  itself.  These 
topics  are  discussed  in  the  following  sections. 

It  must  be  noted  that  the  requirement  to  produce  weather  contours  in  terms  of 
reflectivity  rather  than  power  or  other  units  was  actually  of  minor  importance 
in  this  effort.  This  is  not  to  say  that  the  requirement  is  not  important,  but 
rather  that  the  weather  contour  device  to  be  built  would  be  inherently  accurate 
and  so  not  affect  the  calibration. 

The  calibration  of  the  weather  detection  system  is  overwhelmingly  influenced 
by  the  characteristics  of  the  radar  with  which  it  is  employed.  Effects  of 
sensitivity  time  control  (STC),  circular  polarization  (CP),  and  the  MTI  can 
introduce  errors  considerably  greater  in  magnitude  than  can  be  expected  by 
implementing  any  averaging  or  thresholding  techniques. 
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WEATHER  DISPLAY  ACCURACY. 


One  of  the  measures  of  system  performance  is  the  accuracy  with  which  the 
weather  cell  can  be  displayed.  This  accuracy  can  be  measured  in  two  related 
ways  called  level  accuracy  and  positional  accuracy.  Level  accuracy  is  a 
measure  of  the  variation  in  the  actual  weather  reflectivity  that  causes  a 
contour  to  be  generated  in  a  given  cell.  Positional  accuracy  is  a  measure 
of  the  variation  in  the  generated  contour  from  the  actual  position  of  the 
threshold  crossing  of  the  weather.  The  level  accuracy  is  dependent  on  the 
averaging  done  before  the  threshold  decision.  It  is  not  a  function  of  the 
weather  geometry. 

The  level  accuracy  is  computed  by  considering  the  effects  of  averaging  on  log 
normal  video.  The  standard  deviation  of  the  weather  reflectivity  is  5.57  dB. 
After  averaging,  the  standard  deviation  is  reduced  to  0.36  dB.  This  is  equiv¬ 
alent  to  averaging  over  80  independent  samples,  16  of  which  are  derived  from 
the  range  averaging  and  5  of  which  are  derived  from  azimuth  averaging  (calcu¬ 
lated  in  appendix  A).  The  resolution  cell  size  necessary  to  generate  this 
number  of  independent  samples  is  0.8  nautical  miles  (nmi)  by  1.4°. 

The  positional  accuracy  of  the  system  is  roughly  equivalent  to  the  display 
accuracy.  Unlike  level  accuracy,  the  positional  accuracy  of  the  system  is  a 
function  of  the  storm  characteristics;  namely,  the  gradient  of  the  reflectivity 
(change  of  reflectivity  with  distance).  If  the  storm  has  a  high  gradient,  the 
edges  are  fairly  well  defined,  and  no  problem  results.  If  the  gradient  is  low, 
considerable  error  can  be  introduced  by  small  variations  in  weather  level.  The 
effect  of  a  small  change  in  level  gradient  on  the  presentation  of  a  weather  cell 
with  high  and  low  reflectivity  is  shown  in  figure  1. 

As  more  samples  are  taken,  the  positional  accuracy  decreases.  This  effect  is 
particularly  pronounced  when  the  azimuth  time  constant  is  increased,  because 
the  signals  take  a  long  time  to  decay  and  tend  to  elongate  the  contours  in  the 
azimuth  direction.  This  is  a  common  side  effect  of  integration.  On  observing 
the  display,  the  averaging  parameters  which  seemed  to  produce  the  best  display 
were  0.8  nmi  of  range  averaging,  and  five  independent  samples  in  azimuth. 

DISPLAY  RESOLUTION. 

The  resolution  limits  of  the  display  depend  on  the  operational  requirements  of 
the  air  traffic  controller.  It  is  difficult  to  quantify  this  requirement,  since 
it  depends  not  only  on  the  weather  characteristics,  but  also  on  the  characteris¬ 
tics  of  the  display  and  the  preferences  of  the  controller  using  the  display. 

The  display  resolution  was  chosen  to  produce  a  nearly  square  cell  0.4  nmi  on 
a  side  at  16  nmi.  This  corresponds  to  an  azimuth  dimension  of  about  1.4°, 
which  is  just  slightly  larger  than  an  antenna  beam  width,  which  is  the  minimum 
resolution  attainable  in  azimuth.  Weather  displays  for  radar  typically  use 
dimensions  of  this  order. 
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The  display  cell  size  was  somewhat  smaller  than  the  actual  resolution  of  the 
system.  This  was  because  it  was  felt  that  the  coarseness  of  the  presentation 
using  the  actual  system  resolution  would  seriously  detract  from  the  appearance 
of  the  display.  In  fact,  without  the  deholing  and  character  generation  (to  be 
discussed  later),  the  display  would  probably  be  unusable  even  with  a  resolution 
this  large. 

WEATHER  CONTOURING. 

The  display  of  discrete  weather  information  on  a  radar  display  is  referred  to 
as  contouring.  One  way  to  display  these  data  is  by  thin  isoreflectivity  lines 
in  much  the  same  way  as  altitude  is  displayed  on  a  topographic  map.  This  has 
the  advantage  of  disturbing  very  little  of  the  display  (at  least  in  theory) 
and  contributing  very  little  of  the  total  display  brightness.  In  past  sys¬ 
tems,  the  raggedness  of  the  contours  has  been  the  primary  objection  to  their 
use.  This  raggedness  is  caused  by  the  inherent  fluctuation  of  the  storm  cell 
as  was  discussed  above  in  BACKGROUND.  This  has  largely  been  overcome  in  the 
NAFEC  system  by  improved  averaging  and  processing  techniques.  Another  way 
to  display  weather  is  by  shading  the  interior  of  the  contour  lines.  This 
creates  a  display  in  which  the  reflectivity  of  the  weather  cell  is  indicated 
by  the  brightness  of  the  shading.  This  system  uses  two  levels  of  such  shading 
to  indicate  the  presence  of  weather  above  two  preset  thresholds.  Typically, 
these  thresholds  would  be  chosen  to  indicate  moderate  and  severe  weather. 
Combinations  of  these  two  techniques  are  also  possible  as  in  the  second 
proposal  below. 

ATS  personnel  in  cooperation  with  FAA  Systems  Research  and  Development  Service 
(SRDS)  proposed  the  following  two  contours  for  implementation  in  this  device. 
The  first  contour  type  (called  contour  type  1)  is  the  shaded  version  described 
above;  i.e.,  two  threshold  levels  indicated  by  two  levels  of  intensity  shading. 
The  second  type  (called  contour  type  2)  consists  of  a  low-level  intensity 
shading  to  indicate  the  upper  threshold  level  and  a  continuous  isoreflectivity 
line  to  indicate  the  lower  threshold  level.  An  illustration  of  these  proposals 
is  shown  in  the  first  two  figures  of  figure  2. 

There  are  several  advantages  of  the  first  proposal.  For  example,  there  is  no 
ambiguity  in  deciding  if  any  given  area  contains  weather,  especially  when  the 
contour  line  extends  off  the  visible  area  of  the  PPI,  causing  uncertainty  as 
to  which  side  of  the  contour  contains  weather  and  which  does  not.  The  shaded 
version  is  simple  to  produce,  and  it  ordinarily  would  not  require  as  much 
smoothing  as  the  line  contour,  because  the  eye  tends  to  integrate  jagged  con¬ 
tours  of  this  type  better  than  thin-line  contours. 

The  advantage  of  the  second  proposal  is  primarily  the  absence  of  a  weather 
presentation  over  much  of  the  display  area.  It  was  hoped  that  this  would 
minimize  the  interference  between  the  ATC  presentation  and  the  weather 
presentation. 
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A  third  idea  (contour  type  3)  which  waa  also  implemented  was  the  addition  of  a 
series  of  dots  on  the  weather  side  of  the  proposed  line  contour.  This  technique 
removes  the  ambiguity  when  the  entire  contour  line  was  not  visible. 

SYMBOLIC  REPRESENTATION. 

To  efficiently  develop  a  system  to  generate  a  variety  of  contours,  it  was 
first  necessary  to  find  a  method  for  displaying  the  required  information  as 
a  set  of  smaller,  easy-to-generate  elements.  The  approach  used  is  similar  to 
that  used  to  generate  characters  on  a  cathode-ray  tube  (CRT)  for  alphanumeric 
displays.  Each  range  azimuth  cell  is  decomposed  into  an  eight-by-eight  array 
of  smaller  dots.  This  enables  the  device  to  generate  one  of  a  set  of  special¬ 
ized  symbols  chosen  for  contour  representation.  The  symbols  are  shown  in 
appendix  B  pages  B-l  to  B-9.  This  technique  results  in  the  improved  quality 
of  the  weather  presentation.  To  use  this  technique,  the  contour  device  must 
possess  some  algorithm  to  determine  which  contour  element  to  generate  at  any 
given  time.  This  is  not  a  trivial  problem,  and  no  attempt  was  made  to  solve 
this  problem  analytically.  Instead,  the  "algorithm,"  consisting  mainly  of 
judgment  and  intuition,  was  accomplished  by  hand  for  every  possible  case  and 
then  coded  into  a  read-only  memory  (ROM)  for  the  hardware  implementation.  This 
list  is  shown  in  appendix  B  pages  B-10  to  B-22.  The  principal  advantages  of 
this  system  are  the  almost  complete  flexibility  of  the  contours  possible  and 
the  ease  with  which  they  are  implemented.  Only  64  eight-bit  words  must  be 
coded  to  change  the  contour  type,  and  none  of  these  affect  the  decision  coding. 

DESCRIPTION  OF  EQUIPMENT. 

The  contour  generator  is  built  in  a  19-  by  10-inch  rack-mounted  cabinet  and 
contains  all  of  the  hardware  for  contour  generation.  This  would  normally  be 
separated  into  two  sections,  one  at  the  radar  site  and  the  other  at  the  indi¬ 
cator  site.  The  system  can  accept  triggers,  video,  and  azimuth  information 
from  one  source  and  simultaneously  generate  video  synchronized  to  the  triggers 
and  azimuth  information  from  another,  possibly  different,  source.  The  constants 
of  integration  and  range  averaging  are  adjustable  from  the  front  panel.  Also, 
the  contour  level,  type,  and  the  deholer  function  are  front-panel  controllable. 


THEORY  OF  OPERATION  AND  SYSTEM  TESTS 


This  section  describes  the  functional  operation  of  the  contouring  device.  The 
signal  flow  can  be  traced  on  the  system  block  diagram  shown  in  figure  3.  In 
addition,  a  simplified  schematic  diagram  is  presented  for  each  subunit  in 
figures  4  through  9.  Timing  diagrams  and  other  discript ive  diagrams  are  shown 
in  figures  10  through  15.  Attention  is  given  to  the  processing  performed  by 
each  unit  of  the  contouring  device,  and  no  attempt  is  made  to  provide  a  complete 
technical  description  of  the  type  usually  found  in  instruction  manuals. 
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RANGE  AVERAGING. 


The  range  averaging  process  is  accomplished  in  two  steps.  The  first  is 
integration  over  a  2.4-microsecond  interval  using  an  analog  integrator, 
and  the  second  is  the  digital  summation  of  from  one  to  four  analog  samples. 

The  incoming  signal  is  first  converted  into  a  current  in  the  voltage-to-current 
converter  (see  figure  4).  The  current  derived  in  the  converter  is  used  to 
charge  one  of  two  capacitors,  Cl  or  C2.  These  capacitors  are  charged  alter¬ 
nately  and  sampled  by  the  analog-to-digital  (A/D)  converter  when  they  are  not 
being  charged.  The  timing  for  this  operation  is  illustrated  in  figure  10. 
Notice  that  the  capacitor  voltage  is  stable  during  the  A/D  conversion,  so  a 
sample-and-hold  circuit  is  not  necessary.  The  capacitor  voltage  at  the  end 

of  the  interval  is  given  by  I(t)dt,  and  since  the  current,  I(t),  is 

given  by  I(t)KgV(t),  where  g  is  the  transconductance  of  the  voltage-to-current 
converter,  the  final  capacitor  voltage  is  seen  to  be  VC5S  &.  y(t)dt  the 
integral  of  the  input  voltage  with  respect  to  time.  cjo 

The  problem  with  this  technique  is  that  the  capacitors  are  temperature 
sensitive.  This  effectively  changes  the  gain  of  the  integrator  and  causes 
fluctuation  of  the  generated  contours.  Since  this  only  affected  the  absolute 
calibration  of  the  signal  and  not  the  quality  of  the  display,  these  changes 
were  tolerated.  This  problem  could  be  corrected  by  using  an  improved  method 
of  analog  integration  or  faster  conversion  and  exclusively  digital  averaging. 

The  digital  summation  is  accomplished  by  delaying  and  inverting  the  input 
signal  and  then  summing  this  signal  with  the  input  in  a  recursive  adder  (see 
figure  4).  The  response  of  this  circuit  is  seen  to  be  a  running  sum  of  the 
previous  n  pulses,  where  n  is  determined  by  the  delay  time  introduced  prior  to 
summation.  This  delay  is  varied  by  changing  the  clock  rate  of  the  delay  shift 
register  to  change  the  summation  from  one  through  four  samples  integrated.  The 
timing  for  this  operation  is  shown  in  figure  11.  A  photograph  of  the  processed 
signal  after  range  averaging  is  shown  in  figure  16. 

As  can  be  seen  from  this  photograph,  the  range  averaging  produces  a  smoothing 
effect  only  in  the  radial  direction.  This  can  be  observed  as  a  range  elonga¬ 
tion  of  virtually  all  targets.  Notice  that  there  seems  to  be  very  little 
azimuth  correlation.  This  is  still  not  suitable  for  thresholding. 

AZIMUTH  AVERAGING. 

Averaging  in  azimuth  was  accomplished  after  range  averaging  by  a  conventional 
digital  integrator.  The  schematic  diagram  of  this  device  is  shown  in  figure  5. 
The  integration  cell  size  is  4.8  microseconds,  and  there  are  160  such  cells  in 
range.  For  each  of  these  cells,  the  video  input  is  compared  to  the  current 
contents  of  the  cell.  This  represented  the  exponentially  weighted  average  of 
the  past  cells.  A  fraction  of  the  difference  was  added  to  the  cell  contents. 

The  multiplier  of  the  difference  determined  the  time  constant  of  the  integrator. 


The  multiplier  could  be  one  of  four  values,  1/8,  1/16,  1/32,  and  1/64.  This 
corresponds  to  the  sampling  of  from  15  to  127  effective  samples  as  described 
in  appendix  C.  The  last  two  time  constants  are  for  test  purposes  only  and 
are  not  useful  in  operation,  since  they  correspond  to  averaging  over  four  to 
eight  beamwidths.  The  combined  effect  of  averaging  in  range  and  azimuth 
produces  the  two-dimensional  weighting  function  shown  in  figure  12. 

A  photograph  of  the  system  output  after  azimuth  filtering  is  shown  in  figure  17. 
Notice  that  the  picture  seems  to  be  blurred  due  to  the  low-pass  filtering. 

This  is  evidence  of  data  smoothing  in  the  azimuth  direction.  This  is  much  more 
suitable  for  thresholding  than  the  range-averaged-only  signal.  The  transfer 
function  of  this  filter  is  discussed  in  detail  in  appendix  C  prior  to  the 
calculation  of  effective  samples. 

THRESHOLDING. 

Thresholding  occurred  at  six  independently  set  levels.  These  levels  were 
selectable  to  about  1  dB  on  a  60-dB  dynamic  range  log  scale.  Each  level  had  a 
hysteresis  of  about  1  dB  in  an  attempt  to  reduce  the  noise  contribution  for 
storms  with  very  low  gradients.  If  the  input  signal  is  greater  than  0.5  dB 
above  the  threshold,  a  threshold  crossing  is  declared  by  setting  the  asso¬ 
ciated  threshold  flip-flop.  To  reset  the  threshold  flip-flop,  the  signal  must 
drop  to  a  level  below  -0.5  dB  of  the  threshold  level.  This  is  the  digital 
equivalent  of  the  Schmitt  trigger,  commonly  used  for  the  detection  of  digital 
signals  in  the  presence  of  noise.  Referring  to  figure  6,  the  comparators  set 
the  flip-flop  when  the  signal  input  was  greater  than  the  threshold  setting. 

The  comparators  reset  the  flip-flop  when  the  signal  was  less  than  the  thresh¬ 
old  setting.  When  the  signal  was  equal  to  the  threshold  setting,  no  change 
occurred.  The  threshold  and  input  were  both  quantized  to  five  bits.  It  is 
this  quantization  that  sets  the  threshold  hysteresis.  The  transfer  function 
for  this  circuit  with  a  typical  input  is  shown  in  figure  13.  After  the  thresh¬ 
old  decision  was  made,  the  flip-flop  outputs  were  encoded  to  produce  a  three- 
bit  binary  number  representing  the  highest  threshold  crossed.  It  was  this 
number  that  was  processed  in  the  remaining  sections  of  the  system.  A  photo¬ 
graph  of  the  threshold  output  is  shown  in  figure  18.  The  intensity  indicates 
the  threshold  level.  Notice  that  the  edge  of  the  thresholded  signal  is  rea¬ 
sonably  smooth  and  continuous.  This  would  not  be  the  case  without  prior  range 
and  azimuth  smoothing. 

SCAN  STORAGE  (WEATHER  MAP). 

After  the  signal  is  thresholded,  it  is  stored  in  a  memory  organized  as  a 
1 60-by-256-cel 1  range-azimuth  map.  This  storage  was  necessary  for  realignment 
of  the  averaged  weather  to  superimpose  the  contour  over  the  location  of  the 
weather.  In  practice,  the  weather  contours  appear  slightly  less  than  one 
antenna  scan  late.  This  was  not  a  problem,  because  of  the  slow  moving  nature 
of  weather.  In  the  case  where  some  form  of  slow-speed  data  transmission  is 
used,  the  data  would  show  up  some  minutes  late,  regardless  of  this  effect. 

This  also  provides  an  easy  method  of  aligning  the  video  to  an  unstaggered 
trigger  when  this  is  necessary. 
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The  storage  of  data  in  this  map  depends  on  the  previous  contents  of  the  map 
in  much  the  same  way  as  was  described  in  the  section  on  the  operation  of 
the  azimuth  integrator.  The  difference  is  that,  while  in  the  case  of  the 
azimuth  integrator  the  difference  was  multiplied  by  the  integration  constant, 
in  the  case  of  the  scan  storage,  any  difference  in  level  simply  changed  the 
cell  by  one  unit.  This  provided  a  small  amount  of  scan-to-scan  integration. 

It's  primary  purpose  is  to  reduce  the  disturbance  caused  by  a  momentary  high- 
level  detection  caused  by  a  very  strong  aircraft  target  or  similar  phenomena. 
Additionally,  the  scan  storage  provides  the  means  of  synchronizing  the  input 
and  output  of  the  contour  generator. 

The  input  section  was,  as  far  as  timing  was  concerned,  completely  independent 
of  the  output  section.  In  other  words,  the  input  and  output  sections  could  be 
driven  by  different  azimuth  change  pulses  (ACP’s),  azimuth  reference  pulses 
(ARP's),  and  trigger  pulses.  This,  among  other  things,  provided  the  means  for 
realignment  of  the  displayed  contour  video.  Also,  any  positional  error  caused 
by  processing  delays  could  be  overcome  simply  by  driving  the  contour  generator 
with  the  same  triggers,  etc.,  as  were  used  by  the  display  being  driven. 

In  an  operational  system,  the  data  from  this  memory  would  be  transmitted  via 
phone  lines  or  some  other  means  to  the  indicator  site.  At  this  point,  the 
information  is  read  into  a  duplicate  memory  and  used  to  generate  the  display. 

The  data  are  not  directly  read  into  or  out  of  memory,  but  are  first  loaded  into 
the  memory  controller  along  with  the  desired  address.  This  is  called  a  memory 
request.  In  the  case  of  output  from  memory,  the  request  does  not  contain  data. 
The  output  data  become  available  sometime  before  the  next  request  is  made.  The 
request  services  are  alternated  whenever  both  an  input  and  output  request  are 
made.  The  speed  of  the  memory  is  much  faster  than  the  rate  of  requests,  and 
so  there  is  no  possibility  of  missing  either  input  or  output  requests.  The 
timing  for  this  circuit  is  shown  in  figure  14. 

DEHOLING. 

The  operation  of  deholing  is  a  discrete  integration  process  designed  to  simulate 
the  natural  integration  done  by  the  eye  on  gray,  level  contours.  This  was 
particularly  important  in  the  generation  of  line  contours  which  were  not 
integrated  by  the  eye.  Figure  8  shows  the  simplified  schematic  of  this  device. 
It  is  basically  a  two-dimensional  sliding  window  detector.  The  deholer  received 
the  three-bit  video  from  the  scan  memory  and  immediately  selected  the  contour 
level  to  be  displayed.  This  thresholded  information  was  processed  serially  in 
much  the  same  manner  as  simple  thresholded  video.  This  one-bit  video  was 
delayed  and  sequenced  to  present  the  deholer  area  threshold  section  with  nine 
bits  of  information  representing  a  three-by-three  array  of  weather  cells.  The 
following  rules  were  applied  to  this  array  to  generate  another  serial  video  bit 
stream:  (1)  If  the  center  cell  is  empty,  then  do  nothing;  otherwise,  total  the 
surrounding  cells,  and  (2)  if  the  number  of  surrounding  cells  exceeds  the  preset 
deholer  threshold,  declare  that  the  center  cell  and  all  of  the  eight  surrounding 
cells  contain  weather;  i.e.,  set  them  to  one  (refer  to  figure  15). 
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Notice  that  the  minimum  block  size  is  now  the  three-by-three  array.  This 
had  several  consequences.  First,  all  of  the  gaps  and  holes  in  the  weather 
map  were  filled  in  or  made  smaller.  Second,  all  of  the  isolated  cells  were 
removed.  Third,  some  small  amount  of  information  was  lost  in  the  form  of 
decreased  resolution,  but  this  is  a  natural  consequence  of  the  integration 
process.  The  output  of  the  deholer  is  shown  before  contouring  in  figure  20. 
For  comparison,  figure  19  shows  the  same  data  with  the  deholer  switched  off. 
Notice  that  the  great  majority  of  small  discontinuities  have  disappeared  and 
that  the  edges  which  were  jagged  have  been  smoothed.  This  demonstrates  the 
integration  capability  of  the  deholer  circuitry.  Also  note  that  a  larger  hole 
can  collapse  to  become  a  small  hole.  This  is  an  unavoidable  consequence  of 
the  deholing  operation,  but  should  not  present  any  problems,  because  it  is 
still  an  accurate  presentation  of  the  actual  weather  being  detected. 

CONTOUR  GENERATION. 

It  was  advantageous  to  generate  the  different  required  contours  with  the 
same  equipment  and,  if  at  all  possible,  with  the  same  technique.  This  was 
the  principal  reason  the  symbolic  representation  technique  described  earlier 
was  chosen.  The  implementation  of  this  technique  will  now  be  discussed 
with  reference  to  the  simplified  schematic  in  figure  9. 

The  incoming  video  from  the  deholer  was  sequenced  into  a  three-by-three  array 
in  exactly  the  same  manner  as  in  the  deholer  just  described.  The  array  was 
then  examined  by  the  decision  ROM  to  determine  which  of  the  16  symbols  of  the 
chosen  contour  type  should  be  displayed  in  the  location  at  the  center  of  the 
three-by-three  array.  This  is  accomplished  for  each  range-azimuth  cell  of  the 
display.  This  decision  was  accomplished  independently  of  the  countour  type  at 
this  point  because  of  the  restrictions  in  the  types  of  symbols  that  could  be 
used  to  generate  contours.  This  is  described  in  detail  in  appendix  B.  In 
addition,  this  appendix  also  contains  the  contents  of  the  decision  ROM  with 
the  address  given  as  a  three-by-three  array. 

The  sequence  of  symbols  to  be  generated  was  saved  in  the  symbol  shift  register 
for  the  entire  duration  of  the  display  cell.  This  information  was  fed  to  the 
contour  symbol  generator  which  consisted  of  the  symbol  table  ROM  and  an  output 
shift  register.  The  contents  of  the  symbol  ROM  is  also  shown  in  appendix  B. 

A  total  of  four,  16-symbol  sets  was  stored  in  each  ROM,  and  each  symbol 
consisted  of  an  eight-by-eight  array  of  dots.  The  outputs  of  each  contour 
generator  were  then  sent  to  the  mixer  where  they  were  all  combined  and  out¬ 
putted  to  the  display.  The  final  system  output  for  several  cases  is  shown  in 
figures  21  through  28.  Because  the  symbol  set  can  be  modified  by  changing  the 
contents  of  a  single  ROM  and  because  the  programing  effort  is  almost  trivial, 
this  device  would  make  an  ideal  test  bed  for  the  development  and  test  of  other 
contouring  formats. 

The  line  contours  do  not  become  objectionable  until  the  map  information  is 
introduced.  In  areas  where  map  information  is  not  complicated  or  congestive, 
the  line  contours  may  be  preferable,  but  in  areas  where  the  map  information 
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is  complicated,  the  gray,  shaded  contours  are  to  be  preferred.  Probably 
both  of  these  contours  should  be  offered  to  suit  the  controllers'  preference. 
The  dots  in  contour  type  3  would  be  very  helpful  when  the  weather  extends 
off  the  screen,  because  it  would  better  indicate  the  weather  position.  It 
was  not  observed  that  the  shaded  levels  had  a  detrimental  effect  on  the  detec¬ 
tion  of  weak  radar  targets  in  weather  and,  in  fact,  may  have  enhanced  their 
detectability. 


CONCLUSIONS 


1.  The  subject  weather-contouring  device  is  usable  for  ATC  purposes  and 
meets  the  requirements  specified  for  this  project. 

2.  The  accuracy  of  the  total  weather  detection  system  was  not  limited  by  the 
contour  device  or  the  averaging  that  was  done.  The  most  severe  limitations 
result  from  the  inherent  characteristics  of  the  radar. 

3.  A  measurement  resolution  of  0.8  nmi  and  16  or  32  effective  azimuth  samples 
provided  sufficient  samples  for  an  ASR-type  radar  to  enable  the  contour 
device  to  generate  useful  contours. 

4.  The  display  resolution  of  the  contour  device  was  0.4  nmi  by  16  ACP  (1.4°). 
This  provided  sufficient  detail  for  high-quality  contour  presentation  and,  at 
the  same  time,  provided  sufficient  coarseness  for  proper  action  of  the  deholer 
circuitry. 

5.  The  threshold  of  detection  can  be  set  to  within  1  dB  if  60-dB  log  normal 
video  is  used  as  the  input.  This  accuracy  is  well  within  the  uncertainty  of 
measurement  for  ASR  radars  and  poses  no  threat  of  significantly  inaccurate 
measurements. 

6.  The  use  of  the  deholer  circuitry  significantly  improves  the  quality  of 
the  resultant  display. 

7.  The  contour  generator  effectively  contours  weather,  and  because  it  is 
programmable,  the  display  format  can  be  readily  changed.  The  realignment 
feature  eliminates  the  problems  due  to  azimuth  and  range  shifting,  common  in 
the  earlier  contouring  systems. 


RECOMMENDATIONS 


1.  It  is  recommended  that  the  subject  weather-contouring  device  be  considered 
for  implementation  whenever  a  PPI  display  is  to  be  used. 

2.  It  is  also  recommended  that  this  technique  be  investigated  to  determine  if 
it  has  applications  in  the  all-digital  environment,  including  conversion  of 
the  algorithms  used  to  provide  a  random  write  (vector)  display  output. 
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FIGURE  3.  SYSTEM  BLOCK  DIAGRAM 


FIGURE  A.  RANGE  AVERAGING  SCHEMATIC  DIAGRAMS 


FIGURE  6.  THRESHOLDING  CIRCUIT  SCHEMATIC  DIAGRAM 
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FIGURE  10.  ANALOG  AVERAGING  TIMING 
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FIGURE  II.  DIGITAL  AVERAGING  TIMING 
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THRESHOLD  CIRCUIT  HYSTERESIS 
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FIGURE  15.  DEHOLER  OPERATION 
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FIGURE  19 


RESULT  OF  PROCESSING  SHOWN  BEFORE  DKHOLING 


FIGURE  21.  SYSTEM  OUTPUT,  CONTOUR  TYPE  1,  DEHOLER  ON 
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FIGURE  24.  SYSTEM  OUTPUT,  CONTOUR  TYPE  1,  WITH  MAP 
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APPENDIX  A 


CALCULATIONS  OF  SIGNAL-TO-NOISE  RATIO  IMPROVEMENT  DUE  TO 
AVERAGING  IN  RANGE  AND  AZIMUTH 


The  signal-to-noise  improvement  (defined  as  the  ratio  of  the  variance  of  the 
averaged  signal  to  the  variance  of  the  original  signal)  resulting  from  range 
averaging  a  weather  signal  is,  because  the  range  samples  are  measurements  of 
independent  sample  volumes,  equal  to  the  number  of  samples  averaged.  When 
16  samples  are  averaged,  the  signal-to-noise  ratio  increases  by  a  factor  of  16. 
Since  the  variance  of  log  normal  weather  video  is  31  dB2  without  averaging, 
the  variance  of  the  average  will  be  1.94  dB2,  and  the  root  mean  square  (RMS) 
noise  will  be  1.4  dB. 

The  samples  in  a  given  range  cell  are  partially  correlated  with  time;  that  is, 
they  are  not  independent.  What  decorrelation  exists  results  from  two  major 
sources:  the  scanning  of  the  antenna  beam  and  the  motion  of  raindrops  within 
the  sample  volume. 

The  beam  shape  autocorrelation  function  is  computed  by  assuming  a  Gaussian 
shaped  antenna  beam  (figure  A-l).  The  attenuation  pattern  for  a  Gaussian- 
shaped  beam  is  given  by 

A  =  exp 

where:  A  is  the  two-way  voltage  gain 

0  is  the  angle  off  beam  center 
<p  is  a  constant  related  to  the  beam  width 

expressed  in  dB  this  becomes 

20  log  A  ■  20  log  (expj-02/202  )} 

-  (02/202)  20  log  e 

-  (02/02)  4.34 

when  0  *  1/2  y( y  is  the  two-way  beamwidth)  the  attenuation  is, 
by  definition,  3  dB,  so 

-3  -  (02/02)  4.34  -  (1/2  y)2/02  4.34. 

Solving  for  <f>  in  terms  of  Y  yields 

02  -  .36y2 

0  -  .6y. 

Now,  for  the  case  of  the  ASR-8  radar,  Y  “1  *  or 
So ,  / 

A  *  exp  -(02/2  •  . 36* [l.75  x 

A  “  exp  ~|02(4.56  x  103)|. 


1.75  x  10“2  radians. 
IQ"2]  2} 


A-l 


FIGURE  A-l .  TWO-WAY  ATTENUATION  VERSUS  ANGLE  OFF  BEAM  CENTER 


The  antenna  rotates  at  a  constant  velocity,  w,  so  A  is  a  function  of  time. 
$•*  wt  *  (2rr/rev)(12  rev/min)(l  min/60  sec)t 
0-  1.26t 
6 2«  1.58t2 


so 


A 


exp 

exp 


|-1.58t2  (4.5x103]| 
-|[7.23xl03]t2|. 


The  Fourier  transform  of  this  function  yields 


(n/ 7.23  HO3)1/2  expj-w2/4  x  7.23xlO-5} 
2.08xl0-2  expl-v2  3.46xl0-5}. 


A- 2 


The  power  spectral  density  is  therefore  given  by 
S(w)  ■  4.34x10”^  exp|w2  6.92x10  5j. 

The  3-dB  bandwidth  is  given  by 

1/2  -  exp  -{  w2  6.92xl0“5} 
w2  -  (ln2)/6. 92xl0-5  -  lxlO4 
w  ■  100  ■  2»f 
f  -  100/2 >r  -  15.9  Hz. 


The  beam-shape  autocorrelation  is  thus, 

R(r)  -  1.47xl0-2  exp  |-3.62xl03  . 


The  composite  variance  for  a  weather  echo  is  due  mostly  to  turbulence 

2  ^  ^ 

and  shear  and  was  given  by  Sirmans  and  Doviak  (reference  1)  as  a  *  o*  *  c4 
2  2  v  c  8 

where  a  is  the  total  variance,  o  is  the  variance  due  to  turbulence, 

2  V  C 

and  <?8  is  the  variance  due  to  shear. 

o8  *  (.42  KvR#2  cos(?e]2  m2/sec 
2 

at  m  2  m2/8ec  approx. 


where  Kv  »  4xl0~3  sec"! 

R  -  20  nmi  -  3.22xl04m 

62  »  5* 

de  '  0* 


shear  constant 
range 

two-way  vertical  beamwidth 
elevation  angle 


c?  *  22.3  m2/sec2 
a ^  •  24.3  m2/sec2 


ov  •  5  m/sec. 


crv  ie  transformed  into  frequency  by  the  Doppler  equation  by 


afv 


2  f 
ov  o 


2  x  5  x  2.9x10 
3  x  108 


96  Hz  . 


Assuming  a  Gauss ian-shaped  power  spectral  density,  this  becomes 
Sw(w)  -  e-*2/2  [607] 2 

so  R(T)  -  1.84x10*  e'l1-85*^5]7-2. 

n 

The  autocorrelation  of  the  weather  signal  is 
R^/T)  -  w  e-^l.SSxloS, 
and  the  autocorrelation  of  the  beam  shape  is 
Rss(T2)  "  Be-^.etaiOS. 


The  composite  correlation  is  then 

R88(t)  *  Rw^RggCT)  *  WBe-T^U-S^xlO5] . 

When  this  is  filtered  by  the  azimuth  integrater,  the  output  spectrum  is  given 
by  SQ(w)  ■  Ss(w)Ha(w)H*a(w>. 

The  resulting  correlation  function  is,  withT*  nT 
Roo^nT^  *  Rss(nT>  tha^n^  *  ha(n)]. 

The  bracketed  term  is  given  by 

ha(T)  x  ha(T)  -  I  |  a[l-a]R  a[l-a]n-R 
N  k*o 

■  ^  a2 [ l-a]n  -  a2 [ 1-a ] . 

Using  the  results  of  appendix  C,  the  h(n)  x  h(n)  term  becomes 

h(n)  x  h(n)  ■  4  S  a[l-a]R  a[l-a]n-R 
"  k“o 

“  a2[l-a] . 

A-4 


The  filtered  signals  autocorrelation  is  given  by 
R00(nT)  -  Rss(nT)  x  a2[a-l)n 


2  1  N 

a  -  R00(o)  *  -s  2  R88(kT)  a2[l-a)n-k 

o  N  k-o  n  =  o 

-I?  e-IkT]2  1.89x105  a2{1_a]-k 

"  k-o 

a2  N 

N  e-,l®k^  [l-a]~k. 

When  a  -  -g-  (i.e.,  the  smallest  time  constant  used) 

2.1  1 

a  ■  TT  TT  1 1  +  .95  ♦  .64  +  .30  +  .1  +  .02  .  .  ] 

O  04  N 

2  i  1  i 

^  *  x  x  3  01  s 
°  64  x  6  X  128. 


For  the  case  of  noise  (where  the  autocorrelation  is  1  at  n“o  and  o  elsewhere) 
2  11  1 

n  64  x  6  *  384 

1 

1  . 

Hence,  the  number  of  independent  samples  is  roughly  one-third  the  number  of 

actual  samples.  The  total  number  of  independent  samples  is  thus  -  x  16  x  16 

2  3 

and  so  the  resultant  variance  will  be  decreased  to  -  -  0.36  dB. 
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APPENDIX  B 


CONTOUR  GENERATOR  ROM  PROGRAMING 


A  three-by-three  matrix  represented  by  a  nine-bit  nun&er  is  fed  into  the 
character  decision  ROM. 


Dj 

E 

E 

E 

E 

E 

0 

E 

The  nine-bit  number  represents  the  matrix  as  follows 


8 

7 

6 

5 

4 

3 

2 

1 

0 

The  leftmost  three  bits  represent  the  top  row  of  the  three-by-three  matrix, 
the  center  three  bits  represent  the  center  row,  and  the  rightmost  three  bits 
represent  the  bottom  row. 


If  any  of  the  bits  are  a  logical  "1,"  the  cell  is  "filled,"  showing  that 
weather  is  present.  For  example,  if  the  binary  number  111  010  101  were 
inputted  to  the  character  decision  ROM,  it  would  represent  the  following 
matrix: 


111  010  101 


The  purpose  of  the  character  decision  ROM  is  to  assign  one  of  16  characters 
to  each  of  the  three-by-three  matrices,  of  which  there  are  512  possibilities. 
These  are  all  listed  in  the  following  eight  pages.  Below  each  of  the  matrices 
is  the  three-digit  octal  number  that  represents  it  (OOOg  -  777g). 

Below  the  octal  number  representing  the  matrix  is  a  hexadecimal  number.  This 
number  represents  the  character  that  corresponds  to  weather  in  that  matrix. 
This  number  is  outputted  from  the  character  decision  ROM  into  the  character 
generator  ROM.  The  character  generator  ROM  outputs  an  eight-by-eight  matrix 
which  is  displayed  as  a  character  representing  weather  in  the  original  three- 
by-three  matrix. 
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APPENDIX  C 


EFFECTIVE  NUMBER  OF  SAMPLES  AVERAGED  BY  A  FIRST-ORDER  DIGITAL  FILTER 


FIGURE  C-l, 


BLOCK  DIAGRAM 


From  the  diagram  (figure  C-l) 


vout 

*  z  *vout 

♦  a(Vin 

vout 

-  (1  -  a) 

2  *  vout 

vout 

-  ( l-( 1-a 

)z-M  - 

^out 

-  H(z)  - 

a 

vin 

l-( l-a)z 

'  in 


by  the  inverse  z  transform  h(n)  ■  a[l-a)n,  where  n  is  the  sample  number. 
If  V  -  2anX 


2 

a 

v 


2  v2 
X 

h  x 


*out  ■  2h2(n)o?  ■  o'  Sh(n)2 

in  in 

this  results  in  a  variance  reduction  of 


2  2 
a  „  a  . 

out  .  out 

— = —  ■  J h‘(n),  which,  for  a  uniform  sampling  is  — s — 


N. 


in 


in 


C-l 


Equating  these  two  gives  a  measure  of  effective  samples 


_1 

Me 


for  a 
yields 


■  j  h2(n)  ■  2 
o  o 

_L  .  92 

Ne  1  -  (1-a)2 

_1  .  «2-  _ L 

Ne  l-l+2a-a2  2-a  i”1 


,1  J.  JL  _1 

8*  16*  32*  64 

Ne  ■  15,  31,  63,  127,  respectively. 
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